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ABSTRACT 

In strong gravitational lensing, the multiple images we see correspond to light rays that 
leave the source in slightly different directions. If the source emission is anisotropic, 
the images may differ from conventional lensing predictions (which assume isotropy). 
To identify scales on which source anisotropy may be important, we study the angle S 
between the light rays emerging from the source, for different lensing configurations. 
If the lens has a power law profile M oc , the angle S initially increases with lens 
redshift and then either diverges (for a steep profile 7 < 1), remains constant (for an 
isothermal profile 7 = I), or vanishes (for a shallow profile 7 > 1) as zi — > Zg. The 
scaling with lens mass is roughly 6 oc M^/*^^"'*''. The results for an NFW profile are 
qualitatively similar to those for a shallow power law, with 6 peaking at about half 
the redshift of the source (not half the distance). In practice, beaming could modify 
the statistics of beamed sources lensed by massive clusters: for an opening angle 9 jet, 
there is a probability as high as P 0.02-0.07 (6*301/0. 5°)""'^ that one of the lensed 
images may be missed (for 2 < < 6). Differential absorption within Active Galactic 
Nuclei could modify the flux ratios of AGNs lensed by clusters; a sample of AGNs 
lensed by clusters could provide further constraints on the sizes of absorbing regions. 
Source anisotropy is not likely to be a significant effect in galaxy-scale strong lensing. 

Key words: cosmology: gravitational lensing — quasars: absorption lines — galaxies: 
jets — galaxies: active - gamma rays: bursts 



1 INTRODUCTION 

In gravitational lensing studies, the emission from the back- 
ground source is usually assumed to be isotropic such that 
the appearance of lensed images depends only on the mass 
distribution in the foreground lensing object and the angular 
position of the source with respect to the lens. The fluxes 
of the images can then be written as Fi — Fsrc where 
is the lensing magnification at the position of image i, 
and the source flux Fsrc is assumed to be the same for all 
images. In this case we can interpret observed flux ratios 
simply as lensing magnification ratios: Fj/Fi — fij/ ^i. (The 
one caveat is that if the source is variable, differences in light 
travel times make it necessary to monitor lensed images and 
synchronize the light curves in the source time frame before 
taking the ratio; e.g., Eigenbrod et al. 2005; Kochanek et al. 
2006; Fohlmeister et al. 2007, 2008.) The ability to associate 
flux ratios with magnification ratios underlies many lens- 
ing applications, including using anomalous flux ratios to 
constrain dark matter substructure (e.g., Metcalf & Madau 
2001; Chiba 2002; Dalai & Kochanek 2002; Keeton, Gaudi 
& Fetters 2003, 2005; Chiba et al. 2005), and using lens 
statistics to constrain the mass function and density pro- 
files of galaxies and clusters (e.g., Keeton & Madau 2001; 



Kochanek & White 2001; Takahashi & Chiba 2001; Oguri 
2002; Ma 2003; Kuhlen, Keeton, & Madau 2004; Oguri & 
Keeton 2004; Oguri & Blandford 2009) Q 

However, many astrophysical sources relevant for 
gravitational lensing have some degree of anisotropy in 
their emission. One example is provided by Gamma-Ray 
Bursters (GRBs)0 Numerical simulations (e.g., MacFadyen 
& Woosley 1999) indicate that the local emissivity is a strong 
function of the angle that the line of sight makes with the 
jet axis, and the interpretation of afterglow observations ap- 
pears to confirm this scenario (Perna, Sari & Frail 2003). 
Other interpretations instead suggest that the emission is 
concentrated in a jet with sharp edges and a range of open- 
ing angles that can be as small as a degree scale (Nakar, 
Granot & Guetta 2004). Recent numerical simulations of 
axisymmetric, magnetically driven outflows (Komissarov et 



In lens statistics, the image fluxes are used to determine 
whether multiple images will be detectable and hence whether 
systems will be identified as lenses. 

^ While there are currently no confirmed cases of GRB lensing, 
the possibility has received considerable attention (e.g., Paczyn- 
ski 1986; Mao 1992; Grossman & Nowak 1994; Holz ct al. 1999; 
Nemiroff et al. 2000; Porciani & Madau 2001). 
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al. 2009) have shown that the 7-ray emitting components 
of GRB outflows magnetically accelerated are very narrow, 
with 9jet ^ 1°- Furthermore, the GRB emission is highly rel- 
ativistic, with Lorentz factors F ~ 100-300 (e.g. Piran 2000). 
This implies that each point on the emitting surface is only 
visible to observers within an angle Svicw ^ 1/F ~ 0.2-0.5 
deg. 

Anisotropy in both emission and absorption is also typi- 
cal of Active Galactic Nuclei (AGNs). The fast TeV variabiF 
ity of the blazars Mrk 501 and PKS 2155-304 was interpreted 
by Giannios, Uzdensky & Begelman (2009) as the result of 
compact emitting regions moving with Lorentz factors of 
F ~ 100 embedded within a jet moving at lower speed. Nair, 
Jin & Garrett (2005) proposed a helical jet model to ex- 
plain observations of the gravitationally lensed blazar PKS 
1830—211. Even if the emission is intrinsically isotropic, ab- 
sorption by clouds within the broad line absorption (BAL) 
region can introduce anisotropy into the net flux out of the 
source. Indeed, Chelouche (2003) argued that differential ab- 
sorption along multiple lens sightlines could be detected in 
spectra of lensed BAL quasars, and Green (2006) suggested 
that this effect might explain variability in the broad emis- 
sion lines of only one image of the wide-separation lens SDSS 
J1004-I-4112 (Richards et al. 2004). The presence of smaU- 
scale structure within the AGN outflow is also supported 
by numerical simulations (Proga, Stone & Kallman 2000). 
While the qualitative picture of the BAL region is generally 
accepted, the quantitative details are still rather uncertain. 
Debate over the location of BAL clouds spans some five or- 
ders of magnitude, ~ 0.01-1000 pc (Elvis 2000; de Kool et 
al. 2001; Everett, Konigl & Arav 2002). Very little is known 
about the size of the clouds, although a number of studies 
suggest an upper limit to the size of < 10^^ cm (Baldwin et 
al. 1995; Elvis 2000; also N. Arav 2009, private communi- 
cation). A cloud of this size at a distance of a few parcsecs 
would produce differential absorption on a scale of arcsec- 
onds. 

The key question, then, is whether anisotropy in real 
astrophysical sources is likely to have a significant impact on 
observed strong lensing. In order to answer this question, we 
need to quantify the angle S between light rays as they leave 
the source on their way to becoming the multiple images we 
observe; only anisotropy on a scale < 5 will be relevant for 
lensing. In Sj2]we set up the general problem of lensing of an 
anisotropic source by a spherical mass distribution. In fJ3]we 
study how the angular beam separation 5 depends on the 
lensing configuration for a wide range of lenses with power 
law mass distributions. In ^ we focus on lenses with the 
Navarro, Frenk & White (1996) profile. In !j5]we use Monte 
Carlo simulations to compute the distribution of 5 angles for 
a population of NEW clusters. We discuss the implications 
of our results in ^ We adopt the cosmological parameters 
Ho = 73 Mpc"^ km s-\ D,^ = 0.24, and Qa = 0.76 (Spergel 
et al. 2007). 



2 LIGHT BEAM SEPARATION IN MULTIPLY 
IMAGED SOURCES 

Fig. [T] shows a schematic representation of the gravitational 
lensing geometry with two images, which is generic when 
the source is sufficiently well aligned with a spherical mass 



a, 



ct. 



Figure 1. Schematic representation of lensing with two images. 
A mass profile shallower than isothermal would also have a third 
image that passes near the galaxy on the opposite side from the 
source, but such central images are faint and difficult to detect 
(e.g., Winn et al. 2004) so we focus on the two main, outer images. 



distribution. When the central density profile is shallower 
than isothermal, the image on the opposite side of the lens 
from the source would be accompanied by a third image that 
lies closer to the lens galaxy (Burke 1981), but such images 
are rarely observed (e.g., Winn, Rusin & Kochanek 2004) so 
in our study we refer to the two outermost images. 

As shown in the figure, j3 is the angular position of 
the (unlensed) source on the sky with respect to the optical 
axis (the line connecting the observer and the lens), 6\^2 
are the angular positions of the two main images, and the 
corresponding defiection angles are ai,2. We are interested 
in the angle 5 between the two light rays that emerge from 
the source to produce the observed images. This angle can 
be written as the sum of the angles (5i,2 between the light 
rays and the line of sight from the observer to the (unlensed) 
source. Simple geometry allows us to identify these angles 
as 

5i = |qi| - |6li| +/? , and &2 ^ \&2\ - \e2\ - P , (1) 

so we can express the angular beam separation as 

5= |qi| + |a2| - (l^il + I6I2I) . (2) 

For a spherically symmetric mass distribution, the defiection 
angle is 



(3) 



where M{0) is the projected mass enclosed within angle 9, 
and Di is the angular diameter distance from the observer 
to the lens, and the sign of the defiection angle is chosen to 
match the sign of 9. The source and image positions and the 
defiection angle are linked via the lens equation, 



/3: 



&{9) 



(4) 



where Ds and Di^ are angular diameter distances from the 
observer to the source and from the lens to the source, re- 
spectively. 
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3 POWER LAW MASS DISTRIBUTION 
3.1 The lens model 

To develop a general understanding of how the angular beam 
separation depends on the lensing geometry and the physical 
properties of the lens, we begin with a simple power law mass 
profile. In length units we write M{R) — AK^ with A some 
constant, so in angular units we have M(6) = A{Di9Y . The 
cases 7 = and 7 = 1 correspond to the familiar cases of 
a point mass lens (PM) and a singular isothermal sphere 
(SIS), respectively. The lens equation takes the form 



117-1 



(5) 



where we use the minus sign when 9 > and the plus sign 
when 6 < 0, and the angular Einstein radius is 



f AG A 
\ c2 



(6) 



Note that if we consider the mass within some fixed physical 
radius we have M (x A and hence 6e oc AI^^^-^ 

For 7 < 1 the lens equation ([5]) formally has two 
solutions for all source positions, although when (3 gets large 
the image on the opposite side of the lens is faint. For 7 = 1 
the lens equation has one or two solutions depending on the 
position of the source, while for 1 < 7 < 2 it has one or 
three solutions. For each solution, the corresponding deflec- 
tion angle is given by 



A. 



7-1 



(7) 



3.2 The angular beam separation 

As we study how the angular beam separation depends on 
the lens redshift zi, we want to keep the physical properties 
of the lens fixed, which is why we elected to write the mass 
profile as M — AR'' where j4 is a constant. To facilitate the 
comparison of models with different power law slopes 7, we 
choose the value of A such that the different masses all have 
the same Einstein radius when the lens is halfway between 
the observer and source. 

Fig. [2] shows the angular beam separation as a function 
of lens redshift for a point mass lens (7 = 0), considering two 
values of the source redshift and different values of the angle 
/3 of the source with respect to the optical axis. A striking re- 
sult is the steep increase of 5 with zi . This and other scalings 
can be understood as follows. The t wo images are located 
at the angular positions 61,2 = (/9± •\/4^^~+~/P')/2, and the 



corresponding defiection angles are cii,2 
Equation then yields 

/A 



(A/As) ei/Oi,, 



(8) 



This equation elucidates the trends apparent in the figure. 
First, since 6™ oc D^J^ we see that 5-pM formally diverge^ 
as the lens approaches the source {zi —> Zs and hence Dis 
0). This divergence occurs for all values of the source angle 
/3. Second, it is clear that 5pM — + as the lens approaches the 



^ In practice the small-angle approximation would break down 
before 5 truly diverges. 



observer {zi —> and Dis Ds). Third, when the source and 
lens redshifts and the lens mass are all fixed, 5pm increases 
with /9. Fourth, when the source is well aligned with the lens 
(/3 <^ Oe), the angular beam separation scales with the lens 
mass as SpM oc M^^^ . Finally, in the opposite limit in which 
(3 ^ 9e, 5pM becomes independent of the mass of the lens 
(although this particular case is perhaps less relevant than 
the others because when f5 ^ Be the counter-image is faint). 

Fig. |3] shows the results for other power law profiles, 
including one that is extended but steeper than isothermal 
(7 — 0.5), the isothermal profile (7 = 1), and one shallower 
profile (7 — 1.5). For all cases we fix /3 = 10~* deg, which is 
small enough that we are always looking at situations with 
multiple images. As we shall see, for shallow profiles the 
angular beam separation is not very sensitive to the choice 
of /3. 

Qualitatively, steep profiles (7 < 1) behave in a similar 
way as the point mass case, with 5 increasing monotonically 
with zi and diverging as the lens nears the source. In these 
cases we find that S increases with the source position /3 (the 
case 7 = is shown in Fig. [2] the case 7 — 0.5 is not shown). 

The isothermal profile (7 — 1) represents a transition 
case, with S increasing monotonically with the lens redshift 
and reaching a finite value as — > Zs. This is another case 
we can understand analytically. For (3 < 9e there are two 
images at positions 61^2 ~ 9e i (3. Combining these with 
the defiection angles qi,2 = (A/As)^£; yields the angular 
beam separation 



Ssi: 



201 



/_A 

VAs 



(9) 



An SIS has 6e oc Dis, so when zi Zs we see that S ap- 
proaches a constant. Another general point is that 5 is in- 
dependent of 13, while depending linearly on the mass of the 
lens (5sis on Oe (x M, where we are considering the mass 
within some fixed physical radius; cf. i]3.ip . 

Profiles shallower than isothermal (7 > 1) show a qual- 
itatively different behaviour: as zi increases, the angular 
beam separation 5 initially rises but then reaches a peak 
before turning over and returning to zero as zi —> Zg. While 
there is no simple, general expression for 5 as a function of 
both the power law slope 7 and source position /3, we can 
find an enlightening analytic result for the limit P <^ 9e- In 
this case the images are near the Einstein radius, 81^2 ~ 9e, 
so the two defiections angles are 01,2 ~ A/As^b, and we 
have 5 ^ 29e{Ds/Dis - 1). Then, for As ^ we have 



9e oc -DiV^ and hence 



5 oc D 



1/(2-7) 
Is 



f- 

VAs 



(7_l)/(7_2) 



(10) 



This result clarifies the distinction between steep and shal- 
low profiles: as the lens nears the source, S diverges for all 
7 < 1, approaches a constant for 7=1, and vanishes for all 
7 > 1 . Finally note that in this limit when /3 is small we have 
the following scaling with lens mass: 5 (x 6e (x M^^^'^~''\ 

While we have explicitly shown results for a single value 
of the angular source position (3, we find that the general 
behaviour of 5 with lens redshift and mass holds for other 
values as well. We have demonstrated this analytically for 
the cases of 7 = and 7 = 1. More generally, for steep 
profiles (7 < 1), which have two images for all /3, we find 
that the rise of 5 becomes increasingly steep as (3 becomes 
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Figure 2. The angular beam separation 5 for a point mass lens as a function of the lens rcdshift. Different curves correspond to different 
source positions. The lens mass is M = 3 X IO^^Mq and the source redshift is = 6 (left panel) and Zs = 1 (right panel). The points 
near the bottom of each panel show the Einstein radius as a function of the lens redshift. 




Figure 3. The angular beam separation <5 for a lens with a power law mass distribution M = A(Di6y^ , as a function of the lens redshift. 
Diflferent curves correspond to different power law slopes. The point mass case (7 = 0) has a mass of M = 3 X 10^^ Mq, while in the 
other cases the mass normalization A is fixed so that all lenses have the same Einstein radius when Di = Ds/2. The source redshift is 
Zs = G (left panel) and Zs = 1 (right panel). We fix the source position (3 = 10~^ deg, but note that for 7^1 the dependence of 5 on ,8 
is negligible. 



larger than 6e. For shallow profiles (7 > 1), by contrast, 
there are two images only when (3 is sufficiently small that 
the source lies inside the lens caustic, and it turns out that 
the dependence of 5 on /3 is quite weak. 



4 NAVARRO-FRANK- WHITE MASS 
DISTRIBUTION 

4.1 The lens model 

The NFW profile (Navarro, Frank & White 1996), 



p(r) : 



(r/r.)(l + r/r.)2 ' 



(11) 
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is fully specified by two parameters, written here as a scale 
radius Vs and characteristic density ps- It is customary to 
trade these two parameters for the (virial) mass Mvir and 
concentration c. The transformation goes by way of the virial 
radius rvir, which is related to the virial mass by Mvir = 
(47r/3) Avir pf'vin where Avir is the virial overdensity, and 
the mean matter density is p{z) = 3HQQ,m{l + z)^/8ttG. 
For the virial overdensity we use the fitting formula provided 
by Bullock et al. (2001): Avir ~ (IStt^ + 82a: - 39x'^)/n{z), 
where x = Q,{z) — 1, and Q{z) is the ratio of the mean 
matter density to the critical density at redshift z. Using 
these ingredients, we can finally specify the original NFW 
parameters: Vs = fvir/c, and pa = Sc p, and 



m{x) 



5,{z) 



Avir (a) 



ln(l + c)-c/(l + c) 



(12) 



The two NFW parameters are not entirely independent: 
numerical simulations reveal a relation between the mass 
and concentration of NFW halos. Bullock et al. (2001) de- 
scribe this relation as a log-normal distribution for the con- 
centration whose median value depends on the halo mass 
and redshift. 



od(M,z) = 



l + z 



/ M Y 



(13) 



where M* is the mass of a typical halo collapsing today. For a 
cosmology with as = 0.9 and w = -1, M, = 1.13 x 10^^ M© 
(Kuhlen, Keeton & Madau 2004). The halos in the ACDM 
simulations by Bullock et al. were best described by the 
parameter values Cs = 9.0 and a = —0.13. Although the 
simulations show a scatter, for simplicity we just use the 
median concentration at each mass and redshift. 

The lensing characteristics of an NFW halo are given 
by Bartelmann (1996). The projected surface mass density, 
scaled by the critical density for lensing, is 



k(x) = 2ks 



where x ; 



1 ' 



r/r 



if a; > 1 




(14) 



(15) 



and Ks = psrsS;.jft is a characteristic surface mass density 
in units of the critical density for lensing, 

„2 



(16) 



(17) 

where g{x) = log(2;/2) -I- 1 — f(x). With this notation, we 
can write the lens equation in scaled coordinates as 

4.Ksg{\x\) 



47rG'AAs ' 
The (scaled) mass inside radius x is 

m{x) ~ 2 j dx! x k[x') — Ansgix) , 



y ^ X 



(18) 



where y — P/Oo is the source angle scaled by the angular 
NFW scale radius, ~ rs/D\. An NFW lens has one tan- 
gential and one radial critical curve, whose scaled radii Xt 
and Xr are given by the solutions of the equations 



: 



1 



d m{x) 



dx 



: , 



(19) 



The tangential critical curve maps to the origin in the source 
plane (i.e., it represents the Einstein ring), while the radial 
critical curve maps to the caustic in the source plane, which 
has angular radius j/caus = m{xr)/xr — Xr. If y < ycaus, an 
NFW lens produces three images: one of the main images 
is outside xt, the other is between xt and Xr, and the third 
image (the one we ignore) is inside Xr- All of these scaled 
variables can be converted into angular variables by multi- 
plying by Oo. 

For an image at angular position 6, the defiection angle 

is 



\a{e)\ 



As 



m{\e\/9Q) 



eo 



(20) 



4.2 The angular beam separation 

Fig. |4] shows the angular beam separation 5 as a function of 
lens redshift for certain values of the lens mass and source 
redshift. From a qualitative point of view, the curves re- 
semble the curve for the shallow power law case in Fig. |3l 
This is not surprising, since the NFW profile is shallower 
than isothermal inside the scale radius r^, and the Einstein 
radius is (much) less than the scale radius for all cases of 
interest. 

The 5{zi) curves peak at a redshift which is about half 
that of the source. This is curious but coincidental; it arises 
from the different redshift dependences of the image posi- 
tions and defiection angles. Fig. [S] shows 5 along with its 
four consituents: the image positions 61^2 and defiection an- 
gles di,2. (These are plotted in units of Oq, but that scale 
factor has a weak dependence on zi.) Neither the image po- 
sition curves nor the deflection angle curves peak at ~ Zs/2; 
the image position curves peak at lower redshift, in keeping 
with the rule of thumb that a lens is most effective when it 
is about half the distance to the source. As far as we can 
tell, it is purely coincidental that the different redshift de- 
pendences cause the peak of the 5 curve to be located at 
~ z,/2. 

Another interesting point from Fig.[5]is that the angular 
beam separation is quite insensitive to the source position 
even though its constituents do depend on /3. This is because 
the image position and defiection angle depend on /3 in a 
similar way: as /3 increases, 9i and &i both increase while 62 
and 02 both decrease (in amplitude). The way these terms 
combine to form the angular beam separation means the 
changes largely cancel and leave 5 relatively insensitive to 

To illustrate one way of interpreting the 5 values, let 
us consider lensing of a beamed source. Depending on the 
angular beam separation 5 and the jet opening angifQ ^jot , 
there may be configurations in which we see only a single 
image even though conventional lens theory (which assumes 
an isotropic source) would predict two. For example, in the 
limit that the angular beam separation is larger than the 
jet opening angle {S > 6 jet), there is no way to arrange the 



^ We take Sjet to be the full opening angle, i.e., from one side of 
the jet to the other. 
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Figure 4. The angular beam separation <5 as a function of lens redshift for an NFW lens with three different source redshifts (lejt panel) 
and three different lens masses (right panel). The source position is taken to be /3 = 0.05/3caus, but 5 is not very sensitive to this choice. 
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Figure 5. Image positions (a;i_2), deflections angles (q!i,2) a^nd angular beam separation ((5) for an NFW lens. All quantities are in units 
of Qq (the 'prime' distinguishes them from the corresponding quantities in angular units). The lens mass is Af = 3 X W^Mq, and the 
source redshift is Zs = 6. In the left panel the source is close to the optical axis (/3 = 0.05/3caus), while in the right panel it is close to 
the caustic (/3 = 0.95/3caus). Note that the image positions and deflection angles depend on the source position, but the angular beam 
separation S is almost independent of /3. 



jet such that we could see both images: we could see one 
image or the other, or nothing at all (if the jet does not 
point along either Ught ray), but never both images. In this 
case the probability that we would "miss" one of the images 
predicted by conventional lens theory is unity. 

For S < 6'jot there is some finite probability of miss- 
ing one of the images. Given 5 and 6 jet, we use simple nu- 
merical simulations to consider all possible jet orientations 
and compute the conditional probability that one image is 



missed, given that at least one image is seen. (The condi- 
tional part of the probability just means we do not consider 
cases where the jet is pointing "away" from us so that we 
cannot see anything.) The results are shown in Fig. [6] In the 
limit S ^ Ojat ^ 1 we can derive a useful analytic approxi- 
mation P « (45) /( TT^jct), which is also shown in the figure. 
Combining the probability results with the S values from 
Fig. |4l we deduce that there is some finite probability that a 
beamed source could be "lensed" by a massive NFW cluster 
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Figure 6. For a beamed source with opening angle ^j^t, the 
solid curves show the probability that one of the lensed images 
is missed, given that at least one image is seen, as a function of 
the angular beam separation S. The probability saturates at unity 
for 5 ^ (?jot (see text). The dotted lines show the approximation 
P^(45)/(7rV)- 



in such a way that we miss one of the images, and that this 
probability could be as high as P ~ 0.02-0.07 (^jet/O.S")-^ 
depending on the redshift of the source. This probability 
is largest when zi ~ Zs/2. (By contrast, the corresponding 
probability for a steep power law mass profile would become 
ever larger as — > Zg.) 



5 COSMOLOGICAL POPULATION OF NFW 
LENSES 

The effects of source anisotropy depend on the angular beam 
separation, which in turn depend on the redshift, mass, and 
profile of the lens. For a lens system in which the lensing 
galaxy or cluster is known, it would be natural to study 
source anisotropy using the specific lens configuration. But 
for a system in which the lensing object has not yet been 
identified, it would be important to understand the statis- 
tical distribution of angular beam separations. We now ex- 
amine that distribution using Monte Carlo simulations of a 
population of halos between the observer and source. Since 
we have seen that the angular beam separation is largest for 
massive lenses (which is not surprising), we focus on cluster 
lenses here. We treat the clusters using the NFW profile, 
which of course is oversimplified but serves to give useful 
first estimates. 

We model the cluster population using the mass func- 
tion from Warren et al. (2006), 



n{M,z) = 0.7234 (a 



^^1.625 ^ 0.2538) e 



-1.1982/0 



where p{z) is the mean comoving matter density at red- 
shift z, and aM is the linear density fluctuation on 
mass scale M. In terms of the power spectrum, a|j — 
{l/2n^) J dkk'^ P{k)W{kRf, where P{k) is the matter 
power spectrum, W{kR) is the window function correspond- 
ing the smoothing of the density field (e.g., Peebles 1993), 
and R is the comoving scale corresponding to a mass M — 
(47r/3)pi?^ 

Let us now consider a source at redshift Zs- The proba- 
bility that its light is (multiply) lensed by a mass on its way 
to the observer, is given by the fraction of the sky that is 
covered by lens caustics[f] 



-flcns — . 

where 



f(M,z)dMdz , 



dV , 



(22) 



(23) 



f[M,z) = 7r/3i„, n{M,z) — {z) , 

and dV/dz{z) is the comoving volume at redshift z. We 
can interpret f{M, z) (once it is normalized) as the joint 
probability distribution for the halo mass and redshift, and 
draw from this distribution as follows. First, we obtain the 
redshift distribution by marginalizing over mass: P(ena{^) = 
Al f{M, z) dM, where Ai is a normalization factor. 

We consider the mass range from Mm\n = lO"Af0 to 
Mmax = 3 X IO^^Mq. As we shall see, reducing the mini- 
mum mass would principally affect the distribution of an- 
gular beam separations at the small-5 end, which is not 
so interesting from the standpoint of expected anisotropy 
scales. Our results are not very sensitive to the upper mass 
cut because higher mass halos are exponentially rare. We 
draw a random lens redshift z\ from Pl^^^i^z), and then con- 
struct the conditional probability distribution for the mass: 
Pl^^^{M\zi) — A2f{M,zi), where A2 is again a normaliza- 



. p{z) dlog a J 
M2 dlogM 



(21) 



tion factor. Note that the factor of Tr/^caus in eq- (|23p ensures 
that each halo is weighted by its lensing cross section. 

Once we have drawn a lens mass and redshift, we 
need to draw a random source position. Since we con- 
sider only multiply-imaged sources, we can restrict atten- 
tion to ^ /3 ^ /3caus. We assume a uniform distribu- 
tion of source positions (in 2-d), which is equivalent to 
P{P)d(3 = (2/?//3caus) '^/S- We note that lensing magnifica- 
tion bias can cause the distribution of source positions to be 
non-uniform. However, since the angular beam separation 
for NFW lenses is not very sensitive to /3 (even for /3 on or 
near a caustic; cf. H4.2|l , our results are not very sensitive to 
the distribution of source positions. 

For each set of {zi,M,j3) values we compute the angu- 
lar beam separation 5, and then repeat the process many 
times to obtain the 5 distribution. Fig. [7] shows the results. 
In analyzing the results, we examine one case in which we 
consider all halos above a fixed mass theshold (left panel 
of Fig. [7)1 , and a second case in which we consider all halos 

® If we were considering a particular flux-limited survey, we would 
need to consider not only the fraction of the sky covered by caus- 
tics (the lensing optical depth) but also the fact that lensing con- 
figurations with higher magnifications are easier to detect (mag- 
nification bias). We find, though, that magnification bias does 
not significantly affect our conclusions about the distribution of 
5 values. 
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Figure 7. Cumulative probability distribution for a lensed source at redshift Zj, = 6 to have an angular beam separation greater than 
5. The left panel shows results for different mass cuts, while the right panel shows results for different image separation cuts. 



that produce a lens image separation larger than some value 
(right panel). 

The range of S values typical for the statistical distri- 
bution is somewhat smaller than the values seen in Fig. |4j 
which is not surprising because the cluster population tends 
to be dominated by clusters at lower masses and redshifts 
than were used for the previous example, and both effects 
tend to reduce S. A useful corollary is that, since 5{zi) is 
largely insensitive to Zs (see the left panel of Fig. the 5 
distributions shown in Fig. [7] are not very dependent upon 
the assumed source redshift. 

The typical value of 5 increases with the mass or im- 
age separation threshold, which is also not surprising. If we 
consider wide-separation lenses with A9 > 10", the median 
angular beam separation is 5 « 5", and there is about a 10% 
chance that the angular beam separation would be larger 
than about 10". As discussed below, this range of S values 
corresponds to possible scales of anisotropy in real astro- 
physical sources, which suggests that source anisotropy may 
have non-negligible if not dominant effects in strong lensing. 



6 SUMMARY AND DISCUSSION 

Motivated by the observation that anisotropic emission is 
ubiquitous in astrophysical sources, we have addressed the 
general problem of gravitational lensing of an anisotropic 
source. For simplicity, in this pilot study we have focused on 
spherical mass distributions. 

The slope of the mass profile plays a crucial role in de- 
termining the angular beam separation 5 between the two 
main lensed images, and hence the probability that the ap- 
pearance of the images will be affected by source anisotropy 
(that we might miss one of the images of a beamed source, 
for example). If the mass profile is steeper than isothermal, 
S increases with the redshift of the lens (for a fixed source 
redshift). If the mass profile is shallower than isothermal, S 



peaks when the lens is somewhere in between the observer 
and the source. In the case of an NFW lens, the peak occurs 
for zi ~ Zs/2. In general, the angular beam separation S is 
not necessarily similar to the image separation, A6'. 

By modeling a cosmological distribution of clusters with 
NFW profiles, we have estimated the distribution of angular 
beam separations — and, by implication, the range of source 
anisotropy scales that are most likely to affect lensing. Typ- 
ical values of 5 for cluster lenses lie in the range of a few 
to tens of arcseconds. We remark that while our statistical 
analysis should be instructive, there are some complications 
we have omitted in this pilot study. If baryons steepen the 
central density profile (compared with the pure NFW pro- 
files we have considered; see Puchwein et al. 2005; Rozo et 
al. 2008), that would tend to increase the S values. In general 
we would expect asphericity to change S only by a factor of 
order unity, although that would introduce some new phe- 
nomenology in terms of how source anisotropy would affect 
lenses with more than two bright images corresponding to 
more than two light rays emerging from the source. 

Our results have implications for a range of astrophysi- 
cal sources. Gamma-Ray Bursts are known to be extremely 
relativistic sources, with Lorentz factors F ~ 100-300. This 
implies that each point on the emitting surface is only vis- 
ible to observers within a sub-degree scale. If the GRB 
emitting region is confined within a region < 1° (as sug- 
gested by recent numerical simulations; Komissarov et al. 
2009), GRBs would be excellent candidates to display effects 
of anisotropic lensing. Since GRBs are short-lived sources, 
GRB lensing would not generally result in contemporaneous 
multiple imaging. Rather, if a GRB goes off behind a mas- 
sive cluster, one would naively expect to see a nearly iden- 
tical (modulo some magnification factor) GRB within the 
same instrumental error circle some months or years later. 
We have found, however, that there could be a small but 
finite probability that one of the images could be missed 
because of the angular beam separation between the light 
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rays that form the lensed images. As a specific example, if 
we consider a GRB at Za = 6 lensed by a massive cluster 
at 21 = 1 (i.e., the scenario depicted in Fig. U), and its jet 
opening angle is 6 jet = 1°, there is a 1% probability that 
the second GRB image will not appear. That probability 
doubles if 6 jet = 0.5°, and it increases to about 6% for a 
lens redshift zi = 3. The possibility of missing images ought 
to be incorporated into statistical forecasts of GRB lensing 
(which have heretofore assumed isotropic source emission; 
cf. Porciani & Madau 2001). 

Other sources that have been suggested to be highly rel- 
ativistic are blazars. Giannios et al. (2009) suggested that 
the fast TeV variability observed in two sources can be ex- 
plained as the result of compact emitting regions moving 
towards the observer with Lorentz factors of ~ 100, and em- 
bedded within a jet moving at lower speed. In this scenario, 
the emission from each blob is beamed within a sub-degree 
scale. Depending on the angle at which the blob is moving 
with respect to the line of sight to the observer, there is 
some probability that a lensed blazar might be missing one 
of the images. If that does happen to be the case in lens- 
ing observations (e.g., when a blazar is observed behind a 
large cluster and the number of images appears anomalous) 
it would provide support for this physical picture of blazars. 

Anisotropy in the net flux leaving the source can also re- 
sult from inhomogeneous absorption within the source. This 
is indeed the case for AGNs, in which dense clouds in the 
broad abosorption line region create a highly anisotropic 
absorbtion pattern. The precise location and size of these 
clouds is still a controversial issue. Estimates suggest that 
their size is not larger than about lO^** cm (e.g., Baldwin et 
al. 1995; Elvis 2000), while their location has been placed in 
a range between 0.01 and 1000 pc (e.g., de Kool et al. 2001; 
Everett et al. 2002). These scales fall in a quite interesting 
range for our problem. A cloud of size 10^* cm at a distance 
of 1 pc would subtend an angle of about 7", which is similar 
to the median angular beam separation for NEW lenses with 
image separations Ad > 10". If the cloud were significantly 
closer to the central engine than 1 pc, it would most likely 
cover both of the light rays that correspond to lensed im- 
ages, so the importance of source anisotropy would depend 
on whether there is significant internal structure within BAL 
clouds on scales smaller than 10^* cm. Conversely, if BAL 
clouds are significantly farther than 1 pc, the importance of 
source anisotropy would depend on the covering fraction of 
BAL clouds. 

If an AGN is strongly lensed and there is significant 
differential absorption within the source, that would effec- 
tively cause different lensed images to have different source 
fluxes, which would in turn break the connection between 
observed flux ratios and lensing magnification ratios. If this 
complication is not recognized, it could lead to errors in lens 
models and their interpretation. By contrast, if the differ- 
ential source absorption is recognized, the ability to simul- 
taneously probe multiple lines of sight into the source with 
strong lensing would provide a new way to probe the struc- 
ture of the absorbing medium in AGN, which is still very 
uncertain. This possibility is related to the suggestion by 
Chelouche (2003) and Green (2006) that lensed quasars can 
be used to study small-scale structure in quasar outflows. 
One good would to identify differential absorption would 
be to compare flux measurements at both X-ray and op- 



tical wavelengths (Green 2006). Column densities inferred 
for the absorbers in the broad line regions are A^h ^ lO'^^ 
cm^^, much larger than derived from the UV (Green et al. 
2001; Gallagher et al. 2002). Therefore, one expects that a 
light ray passing through an absorber would have a smaller 
X-ray/optical flux ratio than a non-absorbed ray. 

In summary, anisotropy in sources that are gravita- 
tionally lensed could influence the appearance of the lensed 
images — including whether we even see all the images. The 
effect will be most significant for wide-separation lenses pro- 
duced by cluster-mass objects. Source anisotropy will prob- 
ably not dramatically alter the statistics of GRBs, blazars, 
and AGNs lensed by clusters, but its effect may be non- 
negligible and certainly ought to be considered. If effects 
of source anisotropy can be recognized, they would provide 
a unique opportunity to learn more about the small-scale 
structure of the emitting region of the source. 
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